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Abstract 

We have simulated two-colour four-flavour QCD at non-zero chemical po- 
tential n for quark number. Simulations were performed on 8 4 and 12 3 x 24 
lattices. Clear evidence was seen for the formation of a colourless diquark 
condensate which breaks quark number spontaneously, for fi > /i c ~ 
The transition appears to be second order. We have measured the spectrum 
of scalar and pseudoscalar bosons which shows clear evidence for the expected 
Goldstone boson. Our results are in qualitative agreement with those from 
effective Lagrangians for the potential Goldstone excitations of this theory. 
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I. INTRODUCTION 



QCD at finite quark/baryon-number density at zero and at finite temperature describes 
nuclear matter. Nuclear matter at high temperatures (and possibly densities) was certainly 
present in the early universe. Neutron stars consist of dense cold nuclear matter. RHIC 
and the CERN heavy-ion program promise to produce hot nuclear matter in the laboratory. 
Calculating the properties of high density nuclear matter could predict if and where strange 
matter could be produced. Any method which can be used to determine the properties of 
nuclear matter could be adapted to nuclear physics calculations. 

Finite quark-number density is best achieved by introducing a chemical potential \x 
for quark-number, and using the grand-canonical partition function. Unfortunately, this 
renders the Euclidean-time fermion determinant complex, with a real part which can change 
sign. Standard lattice simulations, which rely on importance sampling, fail in this case. 
Attempting to circumvent these problems by using canonical (fixed quark number) ensembles 
fail except at high temperatures |]J because of sign problems. 

Until a simulation method is found which avoids these difficulties, it is useful to study 
models which exhibit some of the properties of QCD at high fi. Now it is expected that, 
at zero temperature, nuclear matter undergoes a phase transition at fi of order one third 
the mass of the nucleon. It has been proposed that at still higher \i the ground state 
is characterised by a diquark condensate Such a condensate would not only cause 

spontaneous breaking of baryon number, but would also spontaneously break colour. Since 
colour is a gauge symmetry, such breaking is realized in the Higgs mode. Thus nuclear 
matter would become a colour superconductor at high \i. 

For this reason we have simulated 2-colour QCD, i.e. SU(2) Yang-Mills theory with 
fermion matter fields ('quarks') in the fundamental representation of SU(2) co i our , and finite 
fi. As well as having colour confinement this theory does exhibit diquark condensation as we 
shall demonstrate in this paper, but for fi > fi c ~ m n /2, since the diquark 'baryons' in the 
same multiplet as the pions, also have mass m^. For \i < the phenomenon is describable 
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as a rotation of the condensate from the chiral to the diquark direction as predicted by 
effective Lagrangian analyses 0. Unlike in true (3-colour) QCD, the diquark condensates 
are colourless, and the broken symmetry is realized in the Goldstone mode, and there is no 
colour superconductivity, but rather superfluidity, as in liquid 3 He. Despite this, we shall 
later argue that this theory is more similar to 2-flavour QCD than one might think (see 
section 4). 

Since 2-colour QCD has a non-negative determinant and pfaffian, even at non-zero //, 
standard simulation methods can be used. We use the hybrid molelcular-dynamics method 
and simulate the theory with 4 flavours of staggered quarks |§. Pfaffian simulations of a 
4-fermion model including a diquark source term have been reported in ||. We have run 
simulations at a moderately large quark mass and an intermediate gauge coupling on 8 4 
and 12 3 x 24 lattices, i.e. at zero temperature. We measured order parameters including 
the chiral and diquark condensates, the quark-number and energy densities, and the Wil- 
son Line (Polyakov loop). The larger lattice allowed us to observe finite size effects and, 
more importantly to measure the scalar and pseudoscalar meson and diquark masses, which 
include all the potential Goldstone bosons in the theory. Preliminary results from these sim- 
ulations were reported at Lattice'2000, Bangalore ]T0|]. The extension of these calculations 



to finite temperature was reported in a recent letter ||11|| . This work builds on early work 



with 8 quark flavours which presented far less conclusive results | 12| . Previous studies of 
diquark condensation in this model with various numbers of flavours have either used the 
approximation where A = in the updating algorithm [|13J (as does [|12J) ; or been in the 



strong gauge-coupling regime 14 



Section 2 introduces the staggered-fermion lattice port of 2-colour QCD. The results of 
our simulation are presented in section 3. Section 4 gives our conclusions and indicates 
future avenues of research. 
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II. LATTICE 2-COLOUR QCD 



Because in 2-colour QCD fundamental quarks and antiquarks lie in the same rep- 
resentation of SU(2) co i our the flavour symmetry group for Nf flavours is enlarged from 
SUz,(Nf) x SUji(Nf) x Uy(l) to SU(2Nf). The pattern of chiral symmetry breaking is 
SU(2N f ) -> Sp(2N f ) rather than the usual SU L (N f ) x SU R (N f ) -> SU v (N f ) 0,0. 3- 
colour QCD with 1 staggered quark (4 continuum flavours) has a Uy{l) x U(l) e flavour 
symmetry. For 2 colours this is enhanced to U(2) \\12\\ . Chiral or quark-number symme- 
try breaking (the chiral and diquark condensates lie in the same U(2) multiplet) occurs 
according to the pattern U(2) —>■ U(l). 

The quark action for 2-colour QCD with one staggered quark is 

s f = E hmv) + m] X + \\[X T T2X + XT 2X T ]\ (1) 

sites 

where is the normal staggered covariant finite difference operator with \x introduced 

by multiplying the links in the +t direction by e M and those in the —t direction by e _/i . The 
superscript T stands for transposition. Note that we have introduced a gauge-invariant Ma- 
jorana mass A which explicitly breaks quark-number symmetry. Such an explicit symmetry 
breaking term is needed to observe spontaneous symmetry breaking on a finite lattice. We 
shall be interested in the limit A — > 0. Integrating out these fermion fields yields 

Ar 2 A , 

pfaffian ' = Jdet(A*A + A 2 ) (2) 

-A T Ar 2 J 

where 

A=p(fi)+m (3) 

We note that ^^4+ A 2 is positive definite for finite A. Hence the pfaffian never vanishes and 
thus, by continuity arguments never changes sign and can be chosen to be positive. Note 
that det(^4) has been shown to be positive in [[RJ. Denoting the 2x2 matrix in equation |2| 
by M. we have seen that its determinant is the determinant of a positive definite matrix and 
thus can be used directly in our simulations. To do this we define Ai by 



M = 


1 
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r 2 










T 2 







1 



(4) 



so that 



M 



A A 
-A* A 



(5) 



and 



A ] A + A 2 
A4t + A 2 



(6) 



We note that det[A4.^ + A 2 ] = det [^4+^4. + A 2 ]. Because we are now dealing with the ma- 
trix Ai^Ai we can use the hybrid molecular dynamics method with 'noisy' fermions M to 
simulate this theory. Here, although we generate gaussian noise with both upper and lower 
components, we keep only the upper components of the pseudo-fermion field to calculate 
det [A^A + A 2 ], which means that we only need to invert A^A + A 2 at each update. Keeping 
only half the components of the pseudo-fermion field is entirely analogous to keeping only 
the fermion fields on even sites in normal QCD simulations. The square root of equation |2| 
is obtained by inserting a factor of | in front of the fermion term in the stochastic action in 
the standard manner. 

We now give a brief discussion of symmetry breaking in this model. This is covered 
in more detail in [12]. At /i = m = A = 0, the U(2) symmetry will break spontaneously. 



Two directions in which it will choose to break are of particular interest. The first is where 
it breaks to give a non-zero chiral condensate (xx)- There will be 3 Goldstone bosons 
corresponding to the 3 broken generators of U(2). These states and their corresponding 
£7(2) generators are 



0\ X T T2X ~ XT~2X T 



(7) 



^2 X T T~2X + XT2X T 



<73 remains unbroken. The second is where U(2) breaks to give a non-zero diquark condensate 
\{x Tt 2X + X T 2X T )- This time the 3 Goldstone bosons and corresponding generators are 

1 x T ^x + x^x T 

&2 ==> XX (8) 

0"3 X^X - X^2X T 

o\ remains unbroken. When /i ^ only 2 Goldstone bosons remain, x Tr 2 e X + X r 2 e X T an d 
X T T2X — X r 2X T - When in addition m/0, only the latter state remains a Goldstone boson. 
A more detailed study of the patterns of symmetry breaking can be performed in terms of 
an effective Lagrangian for the Goldstone modes as in J7|. The only difference is that the 
effective field, denoted £ in that work, here belongs to a symmetric 2x2 tensor representation 
of U{2) rather than to the antisymmetric 2Nf x 2Nf representation of SU(2Nf) of the 
continuum case. 

One can see the remnant U(l) symmetry when A — > by allowing A to become complex. 
The Majorana mass term in equation [1] then becomes |[Ax T T2X + ^*X T 2X T ]- ^ 2 is replaced 
by |A| 2 in the pfaffian. 

Although the 2-flavour theory would be of more interest in the continuum, we have chosen 
to simulate the 4-flavour theory because this represents a single staggered quark species and 
thus has well defined symmetries and a well defined spectrum at all lattice spacings. Unlike 



the 8-flavour case simulated in [I2J it probably does have a sensible continuum limit. 



III. SIMULATIONS OF 2-COLOUR, 4-FLAVOUR LATTICE QCD 

We have simulated 2-colour QCD with 1 staggered quark species (4 continuum flavours) 
on 8 4 and 12 3 x 24 lattices. The simulations reported here are all at f3 = A/ g 2 = 1.5 which 
is roughly the f3 c for the finite temperature transition on an N t = 4 lattice |I7[]. This first 
set of simulations has been performed with quark mass m — 0.1 in lattice units. These 
simulations are currently being repeated at m — 0.025, where the smaller pion mass will 
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give a richer spectrum of Goldstone and pseudo-Goldstone bosons and where a larger portion 
of the relevant phase diagram should be described by effective Lagrangians. (We have also 
performed some zero temperature simulations at /3 = 1.0 and m = 0.05. This has been 
reported in our finite temperature/finite /x letter |IT[ .) Since we wished to take the limit 
A — > 0, we needed A << m. The values we chose were A = 0.01 and A = 0.02 for m = 0.1. 
(At low /x's we also ran at A = 0.) 

The smaller lattice was used to map out the interesting range of \i values, measuring order 
parameters including the diquark condensate {x Tt iX)i the chiral condensate (xx) — (V"'0)> 
and the number density jo = y 9 g^f 1 ■ in addition to measuring these quantities on the 
larger lattice, we also calculated the spectrum of potential Goldstone bosons. The length 
of each 'run' was 2000 molecular dynamics time units, dt had to be chosen as low as 0.005 
for A = 0.01 and 0.4 < /x < 0.975. In figure [I] we have plotted the diquark condensate as a 
function of /x at each A for both lattice sizes. 

Since we are interested in the limit where the symmetry breaking parameter A —>■ 0, 
we have performed a linear extrapolation of the diquark condensate to A = 0. Note that 
the effective Lagrangian calculations J7J suggest that linear extrapolations are the correct 
approach for A sufficiently small, except at /x c . The results of these extrapolations are plotted 
in figure ^. What we first notice is that for /x < 0.2, the extrapolated diquark condensate is 
small enough that we can believe that it should be zero. For /x > 0.35 it is clearly non-zero. 
The points at /x = 0.3 would appear to show finite size rounding were it not for the fact that 
the 8 4 and 12 3 x 24 points are so close together. We think it more likely that /x = 0.3 is so 
close to the transition that the linear extrapolation has broken down. We therfore conclude 
that the system undergoes a phase transition at /x = /x c ps 0.3 < m n /2 = 0.37622(5). If this 
holds true, the fact that /x c is less than m n /2 would indicate that the diquark 'baryons' do 
not exist as free particles but bind into 'nuclear' matter. Over the range 0.35 < /x < 0.6, the 
condensate increases with a curvature consistent with a critical index f3 qq < 1 (the tree level 
effective Lagrangian analysis predicts the mean field result (3 qq = |). Since the condensate 
starts to decrease soon after /x = 0.6, the scaling region is narrow and it would require 
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FIG. 1. (x Tr 2X) as a function of /j: a) on an 8 4 lattice and b) on a 12 3 x 24 lattice. The arrow 
is at fj, = 
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FIG. 2. (x T T2X) extrapolated to A = as a function of /i on both 8 4 and 12 3 x 24 lattices. 

more points to even try to extract this critical index. For fi > m n , the condensate starts 
to decrease, approaching zero for large fi, which would appear to be a saturation effect. In 
fact figure § suggests that the condensate vanishes for \i > n a w 1, indicating a saturation 
phase transition. We shall have more to say about this later. Note that this decrease in 
the condensate with /i is not predicted by the effective Lagrangian analysis, which is only 
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expected to be valid for small /i and m w . This is not surprising if it is indeed a saturation 
effect. Saturation is a result of the fermi statistics of the quarks, and should not be seen in 
a model which only considers the system's bosonic excitations. 

We now turn to a consideration of the chiral condensate, (xx) = (V''0) • This is plotted 
in figure ^ for our 2 different lattice sizes. In the A — > limit, this is expected to be 
constant at its /i = value for /i < /x c . These plots are consistent with this expectation. 
Above /i c , effective Lagrangian studies predict that the condensate merely rotates from 
the chiral direction to the diquark direction, so that the magnitude of the condensate, 
i.e. J (xx) 2 + (x Tr 2X) 2 ; should remain constant and independent of A. Since the diquark 
condensate increases up to fi ~ m f , this means that the chiral condensate should fall, as it 
does. It, however, continues to fall past this point, because of saturation effects, appearing 
to vanish for fi > fi s . To test the prediction that the principal effect for /i < is merely 
a rotation, we have tabulated yj (xx) 2 + (x Tt 2X) 2 i n table |. We see that, for /x < 0.6, the 
magnitude of the condensate is approximately constant and independent of A, so that the 
main effect is a rotation of the condensate from the direction of chiral symmetry breaking 
to that of quark-number breaking. For /i > 0.8, it decreases with increasing \x approaching 
zero for large //. 

In figure [| we plot the quark- number density jo = ^ 91 ^ S/ as a function of /i for our 2 
lattices, jo is consistent with zero as A — > 0, for \i < \i c . (Note effective Lagrangians predict 
that it vanishes quadratically with A in this region.) Above /i c it increases with increasing 
fi approaching the saturation value of 2 (1 staggered quark field of each colour/site, the 
maximum value allowed by fermi statistics) for large /i. Note that figure £| suggests that 
jo = 2 for fi > /i s with fj, s consistent with that predicted from the condensates. This indicates 
that it is this saturation which causes the condensates to vanish for fi > fi s . It also tells us 
that this transition is an artifact of the finite lattice spacing and would recede to infinite \x 
in the continuum limit. 

On the 12 3 x 24 lattice we have measured the propagators for all the potential Goldstone 
bosons, both scalar and pseudoscalar using a single noisy point source for the connected 
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FIG. 3. (ipip) as a function of fi: a) on an 8 4 lattice and b) on a 12 3 x 24 lattice. 
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TABLE I. Magnitude of the condensate y (xx) 2 + (x Tt 2X) 2 as functions of fi and A on a 
12 3 x 24 lattice. Errors are not quoted but they are in the next or a subsequent digit after the 
least significant quoted digit. 



propagators and multiple (5) noisy point sources for the disconnected propagators. We have 
measured both diagonal and off-diagonal zero momentum propagators. 

In figure || we show the masses from fitting the propagator for the diquark state produced 
by applying the operator % T r 2 x — X T 2X T to the vacuum, to the form 

P G {T) = A {exp[-m G T] + ex V [-m G (N t — T)]} , (9) 

valid for large temporal separations T. As argued in section 2, this is the expected Goldstone 
boson for /i > fi c and A = 0. The effective Lagrangian approach indicates that, in the low fi 
phase, this mass should be quadratic in A, while in the high fx phase, it should vanish as \/A. 

12 



SU(2) N f =4 |8=1.5 m=0.1 8 4 lattice 



2.0 



1.5 



1.0 - 



0.5 



o 
x 
o 



A=0.00 
A=0.01 
A=0.02 



e 



bjiL 



_i i i i i i i i i i i i i i i i_ 



o o 1 » 1 » 
0.00 0.25 0.50 0.75 1.00 1.25 

(a) M 

SU(2) N f =4 0=1.5 m=0.1 12 3 x24 lattice 



2.0 



1.5 



1.0 



0.5 



o 

X 

o 



A=0.00 
A=0.01 J 
A=0.02 J* 



8 



_i i i i i i i_ 



0.0^ 

0.00 0.25 0.50 0.75 1.00 1.25 

(b) M 
FIG. 4. Quark-number density as a function of fi: a) on an 8 4 lattice and b) on a 12 3 x 24 



lattice. 



13 



SU(2) N f =4 0=1.5 m=0.1 



GO 
W 

c6 




0.00 0.25 0.50 0.75 1.00 1.25 

/a 

FIG. 5. Mass of the diquark state expected to become a Goldstone boson in the diquark 
condensed phase, as a function of \x for A = 0.1 and A = 0.2. The line is the expected linear 
behaviour expected for the symmetric phase. Also included are the point at A = 0, fj, = and 
A —* extrapolations. 



We have performed extrapolations based on this and plotted them in the figure. For the 3 
points closest to the transition \i = 0.3, 0.35, 0.4, we have performed both extrapolations. 
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For ji < m n /2 and A = 0, we expect 

mc = rn^ — 2/x (10) 

with = m n (fi = 0). We see that for fi < 0.3, the points obtained from quadratic 
extrapolation lie on this straight line, the point at [i = 0.35 is close to the line and that 
for fi = 0.4 lies above this line. All points for 0.35 < // < 1.0 obtained by square root 
extrapolation are negative. The points at = 0.35, 0.4 show the most significant negative 
values (although both are greater than —0.05), and should be considered as transitional. The 
other points are so close to zero that one can easily attribute the difference as being due to a 
combination of the systematic errors in extracting masses from point source propagators and 
higher order terms in the extrapolation. Thus our results are consistent with a transition to 
a phase with a massless (Goldstone) scalar diquark at /i ~ = 0.37622(5), as expected. 

Finally the large value of this diquark mass at fi = 1.2 is further indication that the system 
has passed through the saturation phase transition. 

We now turn to consideration of the other potential Goldstone bosons of the theory. 
Since at this quark mass, is relatively large, we concentrate on those states which are 
expected to have masses of order m n or less over the range of interest. For this reason, 
we consider the pion itself. For /i < /i c , its mass should remain constant at m n , and it 
is created by the operator x e X- F° r A* > /■*« it is expected to mix with the pseudoscalar 
diquark state created from the vacuum by the operator x r T 2 ex + X T 2^X r ■ Rather than trying 
to diagonalize the propagator over these 2 states, we have instead calculated the diagonal 
propagators for each of these states separately. These we fit to the form 
P n {T) = A {expf-m^T] + exp[-m w (N t - T)]} + B(-l) T {exp[-m 61 T] + exp[-m 6 i(iV t - T)}} 

(11) 

or the form with B = 0, and a similar form for the diquark state. These masses are plotted 
in figure |6|. Note that these 2 graphs are almost identical. For fi > fi c , this is expected, since 
the 2 states mix. For /i < /i c this is not, a priori, expected since the diquark state should 
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FIG. 6. Pion (a) and pseudoscalar diquark (b) masses as functions of ji. 
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exhibit a linear decrease with fx for A = 0. However, since A 7^ this gives a small mixing 
with the pion state, which is then the lowest mass particle contributing to this pseudoscalar 
diquark propagator in this region. It is precisely because this mixing is so small that the 
errors are so large. The final indicator that this is the correct interpretation was that the 
fx = 0, A = propagator yielded no sensible mass fit, a sure indicator that the mass was 
large where there can be no mixing. The pion mass fits in this low fx domain are consistent 
with the expectation that the A = pion mass is independent of fx in this region. For fx > fx c 
these masses fall faster than the expected m^/2/x. At least part of the reason for this more 
rapid falloff is the saturation effect. Note that the reason that this mass tends to zero at 
large fx (before the saturation transition) is because the pseudoscalar diquark would be a 
Goldstone boson if m = 0, and as fx increases the relative importance of m diminishes. The 
fact that the ir mass becomes more poorly determined for large fi is because here, the lowest 
mass state in this channel is predominantly the pseudoscalar diquark. 

It is interesting to contrast our findings with those of a study of adjoint quarks in 2-colour 
QCD with A set to zero, implying strictly no mixing; in this case was found to rise 
as « 2fi for fi > fi c , with the signal becoming appreciably noisier in the high density 
phase. Both behaviours are in accord with the predictions of chiral perturbation theory for 
a meson formed from quarks with a symmetric combination of quantum numbers under the 
residual global symmetry, the difference arising due to the distinct Dyson indices of each 
model f7j. 

IV. CONCLUSIONS 

Two-colour lattice QCD with one staggered quark species (4 continuum flavours) has 
been studied at finite chemical potential fx for quark number. We have shown conclusive 
evidence that this theory undergoes a phase transition to a phase characterized by a diquark 
condensate which spontaneously breaks quark-number, at fx — fx c ~ m n /2. This transition 
appears to be second order. The simulations were performed at an intermediate value of the 
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coupling (close to f3 c for the finite temperature transition at /x = on a lattice with N t = A). 
Because of this, the relevant symmetry was the lattice flavour symmetry U(2) rather than 
the SU (8) of the continuum 4-flavour theory. We have presented convincing evidence that 
the scalar diquark is the Goldstone boson associated with the spontaneous breaking of quark 
number for fj, > fi c . In addition we have measured the lightest mass in the pion channel as a 
function of /i. For /i < the chiral condensate and diquark condensate are well described 
as a single condensate of constant magnitude which rotates from the chiral direction for 
/j, < fi c towards the diquark direction as \x is increased above /i c . Despite the fact that the 
quark mass m — 0.1 was large enough that m n was not small relative to the 2-colour QCD 
scale, we saw good qualitative agreement with previous calculations in terms of effective 
Lagrangians (of the chiral perturbation theory form) for \x ^$ m n . This includes the fact 
that the quark number density increases from zero at \i = fi c . 

As fi becomes large, this relationship with effective Lagrangians is no longer valid. Both 
chiral and diquark condensates decrease towards zero as [i becomes large. However, because 
at least part of this is due to the fact that jo saturates at 2 fermions/site, a finite lattice 
spacing artifact, and because we see large finite size effects in the condensates for these /i 
values, studies at smaller lattice spacings as well as on larger lattices would be needed to 
determine how much of this observed high /x behaviour is real. 

We are currently extending these calculations to lower quark mass (m = 0.025) where 
we should be able to measure the complete spectrum of Goldstone and pseudo-Goldstone 
bosons. For this m, the pion mass should be half that at m = 0.1 and the assumptions of 
the effective Lagrangian approach should have more validity, allowing a more quantitative 
analysis. In addition, there is a larger range of n/m^ between [l c and the turnover point. 
This means we can hope to measure the critical index f} qq for this transition. (Preliminary 
attempts to extract this index from short runs at relatively strong coupling were reported 
in our letter on finite T and fi.) 

Let us now compare 2-colour QCD with 3-colour QCD, at finite \i. Most of this discus- 
sion is condensed from the published literature [f§-fj. Since for 2-colour QCD, the diquark 
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condensate is a colour singlet, the spontaneous breaking of quark number is realized in the 
Goldstone mode. This contrasts with true (3-colour) QCD where the condensate is, of neces- 
sity, coloured and the symmetry breaking is realized in the Higgs mode. To see similarities, 
we consider the mode of symmetry-breaking for normal QCD with 2 light quark flavours. 
Here the expected condensate is a flavour singlet and a colour anti-triplet, i.e. antisymmetric 
in both flavour and colour. The pattern of colour breaking is SU (3) x U(l) q — > SU (2) x U(1)q 
where q and Q refer to quark number before and after the breaking. The 5 gluons corre- 
sponding to symmetries broken by the condensate, gain masses via the Higgs mechanism 
by combining with the 5 would-be Goldstone bosons associated with colour /quark-number 
breaking. With respect to the remnant 577(2) colour symmetry, the condensate is a flavour 
and colour singlet. This is precisely the condensate that would be formed from the 2 quarks 
which are in a colour doublet with respect to this unbroken SU(2) (the third quark is a 
singlet with respect to this group) if we ignored the interactions of the gluons which gain 
masses due to the Higgs mechanism. Since the gluon masses produced from the Goldstone 
bosons produced in this manner are of order ^^(Aqcd)^, ignoring such light particle in- 
teractions is presumably a rather drastic approximation. However, one might hope that it 
has at least qualitative validity. (Note that /„. is of the order of 100 MeV, a s (kQ C D) ~ 1, 
giving gluon masses of the same order of magnitude as previous estimates.) Since this rem- 
nant SU (2) theory is also at finite \i this condensate spontaneously breaks only 1 symmetry, 
quark number. As seen above, the Goldstone boson associated with this breaking gives mass 
to one of the gluons via the Higgs mechanism. 

Finally, we note that with an appropriate reinterpretation of the fields, 2-colour QCD 
with a finite chemical potential for quark number can be reinterpreted as 2-colour QCD with 
a finite chemical potential for isospin. It is clear that one can simulate true 3-colour QCD 
for 2 flavours with a finite chemical potential for isospin. Such a program is underway. The 
2-colour theory can be used as a guide to the pattern of symmetry breaking and what to 
expect. 
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